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The optimal entanglement manipulation for a single copy 
of mixed states of two qubits is to transform it to a Bell di- 
agonal state. In this paper we derive an explicit form of the 
local operation that can realize such a transformation. The 
result obtained is universal for arbitrary entangled two-qubit 
states and it discloses that the corresponding local filter is 
not unique for density matrices with rank n = 2 and can be 
exclusively determined for that with n = 3 and 4. As illustra- 
tions, a four-parameters family of mixed states are explored, 
the local filter as well as the transformation probability are 
given explicitly, which verify the validity of the general result. 

Entanglement manipulation is an important issue in 
the studies of quantum information theory. On the one 
hand, it is related to the basic problem of which tasks one 
can accomplish with a given resource of entanglement JjJ- 
. On the other hand, since most applications of quan- 
tum information theory require the maximally entangled 
state for faithful transmission of quantum data, it is nec- 
essary to develop the special technique of entanglement 
manipulation which uses local quantum operations and 
classical communication (LOCC) to produce states with 
the amount of entanglement as great as possible from 
partly entangled states ^,|) . 

There has been much attention recently concerning the 
entanglement manipulation of a single copy of mixed two- 
qubit states ||- p4| . A basic measure of entanglement of 
mixed states is the entanglement of formation, which is 
intended to quantify the entanglement resources needed 
to create the states (JT5). Given a two-qubit density ma- 
trix p, the entanglement of formation is defined as 



E(p) 



(1) 



where the minimum is taken over all decompositions of 
p into pure states 



P = ^Pi\^i){^i\ 



(2) 



and E(tpi) is the entropy of entanglement of the pure 
state \ipi}- Although the expected entanglement cannot 
be increased by any LOCC actions, one can gamble, with 
risk of failure for some probabilities, on obtaining states 
with more entanglement from less entangled states [Q. 
It was proven by Kent etc. jllj that under LOCC oper- 
ations, the entanglement of formation of a single copy of 



two-qubit states exists an upper bound, the maximum 
extractable entanglement, and the corresponding state is 
a Bell diagonal state which is unique up to local unitary 
transformations. In this stage, a critical problem is to 
calculate the detailed LOCC operations and determine 
the probability of success of such a transformation for 
given two-qubit states. 

The local operation on the two-qubit state can be stud- 
ied in a real matrix parameterization representation of 
the state, i.e., the Horodecki representation full . In a 
recent paper [T^ ] , it was shown that local operators per- 
formed on two qubits correspond to left and right mul- 
tiplication by a Lorentz matrix, respectively, and useful 
results have been derived from such a point of view. How- 
ever, until now, the above problem is still open and the 
distinct form of the LOCC actions that realize the op- 
timal transformation is lack. In the presented paper we 
investigate this problem in a particular representation, 
the "Wootters representation" of two-qubit states |l7| . 
We shall show, by introducing an associated operator — 
a squared quantity of the local filter, the explicit form 
of the optimal transformation as well as the probability 
can be derived for arbitrary entangled two-qubit states. 
Moreover, the result displays that the corresponding lo- 
cal filter is not unique for density matrices with rank 
n = 2 and can be exclusively determined for that with 
n = 3 and 4. As illustrations, a four-parameters family 
of mixed states shall be explored, the local filter and the 
transformation probability are given explicitly, so that 
the validity of the general result is verified. 

Main result. — Consider the local operations on a single 
copy of two-qubit states. Recall that the local operator 
can be written as a unitary transformation multiplying a 
Hermitian filtering operator 10 1. One can, without loss 
of generality, describe the transformation as 



P 



fA ® fBpfA g ft 

Hf% ® Pbp) 



(3) 



Here, we have neglected the local unitary transformation 
since they shall preserve the ent ang lement. The local 
filters /a and fg are given by |Io| , |l4 



Sa = 



1 



-(7 2 + am • a), 



1 + a 

fB = -r—rih + bn-a), 
l+o 



(4) 
(5) 
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where m, n are unit vectors and the parameters a and 
b satisfy < a < 1, < 6 < 1. Above expressions have 
excluded the trivial cases of the niters with the deter- 
minants det(/2 — /1(b)) > which differentiate from the 



f 



operators (ji)) and (j|) only by a coefficient less than 1. In- 
troduce an associated operator of the local filter Ja <B> /b 



Fab 



fl 



fl 



det f A det f B 



(6) 



This quadratic form Fab is equivalent to the local filter 
/a<S>/b m a sense that once is provided, the filtration 
parameters a, b, m, and n can be uniquely determined 
from Eqs. (@)-(@)- It can be easily checked, from the 
relation f a(b)I A(B) = det f A (B)h, that the operator F AB 
has the property 



Fab Fab — h ® h, 



(7) 



where the tilde "~" of an operator is defined as Fab = 
(72 <8> (T 2 f ^4^5 c 2 <8> 02. In the following, we shall determine 
the detailed form of Fab for the optimal entanglement 
manipulation which transforms the two-qubit state p into 
a Bell diagonal state. 

We make use of the so-called "Wootters' representa- 
tion", i.e., the particular decomposition of the two-qubit 
state lUl 



random local density matrices tr ap' — trBp' = 5^2, thus 
is Bell diagonal. To prove the converse statement, we use 
the fact that the Bell diagonal state has the property pj] 



trR(p)=J2 x i = L 
Now suppose the transformed state 

, _ E"=l fA ® fB\Xi)(Xi\f A <S> } B 



(13) 



P 



Hpfl ® /I) 



1=1 



(14) 



is Bell diagonal. Here, is given by Eq. (|llj) andp; is 
the probability of it in the decomposition. Direct obser- 
vation can be found that the states \Ei) (i = 1, • • • , n) 
have tilde inner products (Ei\Ej) = (Ei\Ei)Sij with 
(Ei\Ei) positive real numbers. Thus we have 



Y,trR{p')=Y,Pi( E i\Ei) = l- 



(15) 



P=^\X 1 ){2 
1=1 



(8) 



Here, n < 4 is the rank of p. The states \xi) satisfy 

(Xi\xj) = XiSij, \Zi) =0-2(81 cr 2 |x*), (9) 

and the AjS are the eigenvalues, in decreasing order, of the 
Hermitian matrix FL{p) = -\J \fpp~\fp- A combination of 
these symbols defines the concurrence for the entangled 
state p: C(p) = \\ — A 2 — A3 — A4. It turns out that 
the above representation is very useful to determine the 
operator Fab- Namely, we have 

Proposition. — An entangled two-qubit state p given by 
(||) with A„ > can be transformed to a Bell diagonal 
state by the local filter f A <S> /b if and only if the associ- 
ated operator Fab satisfies 

{xi\F AB \xi) = (Xi\xi), i = l, (10) 

Relation ( |Io|) implies that the normalized states 

fA ® fB\Xi) 



\Ei 



\/{Xi\f A ® f B \Xi 



i = l, 



(11) 



have the "tilde inner products" (Ei\Ei) = 1, thus \Ei) = 
\Ei). This yields explicitly 



(x l \F AB \xj) = (xi\xj) =0, i ^ j. 



(12) 



The sufficiency of the proposition can be shown di- 
rectly. Following the fact that the states \Ei) (i = 
1, • • • , n) have the concurrence 1, that is, they are a set of 
maximally entangled states of qubits A and B, one can 
conclude that the transformed state p' has completely 



This yields that (Ei\Ei) = 1 for all i = 1, • • • , n, then the 
relation jic|). 

The analysis presented above also indicates that the 
two-qubit states with A„ = cannot be transformed into 
Bell diagonal states by any local actions efficiently. Note 
that A„ = can only be possible when n — 3 and n = 2 
fLSfl . In such cases, \x n ) is a factorizable state of qubits 
A and B, and cannot be transformed to the maximally 
entangled state (corresponding to \E n ) in the above de- 
scription) by any local filtering operations. The best one 
can do, is to reduce the probability of this term to the in- 
finitesimal. This process, as has been pointed out in Ref. 
|]l3| , transforms the state asymptotically to a Bell diag- 
onal state with a lower rank and an infinitesimal proba- 
bility. 

Equations (|TJJ) and ( |l2] ) suggest a distinct way to con- 
struct the associated local operator F AB . Knowing that 
the rank n = 1 corresponds to pure states and they can 
be purified by the Procrustean protocol Q , we discuss in 
the following the cases of n = 4, 3, and 2, respectively. 

1. n = 4. In this case, the states {\xi),i = 1, • • • ,4} 
compose a complete set of non-orthogonal bases of the 
4-dimensional Hilbert space of the two qubits. Equa- 
tions ( [Toj ) and ( |l2| ) thus provide a clear representa- 
tion of the operator F AB in these bases. Recalling 
{\xi),i = l,---,4} the set of biorthogonal bases to 
i = 1, • • • , 4}, we can present explicitly 



Fab = Y, 



(16) 



To derive detailed expressions for the local filters f A and 
f B , we need to factorize the operator given above. This 
shall be shown later in this section. 
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2. n = 3. Apparently, there exists a unique state \xi) 
in the 4-dimensional Hilbert space which is orthogonal to 
the three linearly-independent states {\xi),i = 1,2,3}. 
Noticing the relation \Ei) — \Ei) (i = 1,2,3), one can 
find that 



The states {|l/i)} have orthogonal normalized tilde inner 
products 



(vi\yj) = s ij> *>i = !>•••, 4. 



(24) 



Construct a new decomposition of the operator Fab — 



{x 4 \F A B\xi)={x i \F A B\xi)=Q, z = 1,2,3. (17) Ei=i l«i><«*l with 

Suppose {x4\Fab\xh) — t/&. The operator Fab can be 
expressed as 



\Xj)(Xj\ 
(Xi \Xi ) 



t 4 4 



\X4) (X4 I 
(xAxi) 2 



(18) 



Substituting this expression into Eq. (Q), one obtains 
t44 = (X4IX4). Thus in this case the local operator Fab 
is also unique and can be expressed uniformly as ([l6]) . 

3. n — 2. Similarly, provided the states 
l^)}; one can easily construct the complete set {\xi),i — 
1 , • • • , 4} in the 4-dimensional Hilbert space of the two 
qubits. In this case, the states \xq) and (0:4) can be cho- 
sen variously. Suppose Uxi), i — 1, ■ • ■ , 4} is one of the 
selected sets satisfying (KJ). Based on an analysis simi- 
lar to the previous case, we can present the form of the 
Hermitian operator Fab as 



F 



AB 



E 

i=l 



(x{ I Xij 



i=3 



\xj)(Xj\ 
\X% \Xi) 



\x 3 ) {Xi I 



\X4){X3,\ 



-). (19) 



{x 3 \x 3 ){x 4 \x 4 ) {X3,\X 3 ){X4\X4) 

Direct calculations from the restriction (^) can show that 

^34 = 0, t 33 = (x 3 \x 3 ), < 44 = (a; 4 |i 4 ). (20) 

This yields again the form of (|l6|). The operator Fab 
now is not unique due to different selections of the states 
{^3), 1^4)}- Specifically, they can be obtained through 
constructing the states as 



ci\x 3 ) 



di\xi) 



C2\X 3 ) 



+ 



d 2 \x 4 ) 



with the complex parameters Ci and di satisfying 
c*c* 2 +dld* 2 =0, 



(21) 



*\2 



1. 



( C ;) 2 + (^)^ = i 



*\2 



(22) 



We now need to show the operator Fab presented in 
Eq. ( |l6| ) can always be factorized as a product form 
of local operators of qubits A and B. To this end, we 
subnormalize the states {|xi), i — 1, ■ • • , 4} such that 



Fab =^2\Vi){Vi\, \Vi. 



\J (Xi\Xi) 



(23) 





-b 


:\yi) 


+ - 


-Kb) 


+ I1/4))], 


*2) 




:\yi) 


+ Iw)) ■ 


-i(lwa) 


+ I1/4))], 


zs) 


-\* 




- Iw)) - 


"I" *(|W3> 


- |»4))], 




"It 


(\yi) 


- Iw)) 


-« 


' - I1/4))]. 



(25) 



One can check directly that the four states \zi) (i — 
1, • • • , 4) have the concurrence zero, thus are factorizable 
states of qubits A and B. In detail, there are 



(zi\zj) = 5 2j , (z 2 \zj) =Sij, 
(z 3 \zj) = S 4j , (z 4 \zj) = S 3j . 



(26) 



Now we show that these relations together with the 
linearly-independent property of the states {\z.i),i = 
1, • • • , 4} can predict that the operator Fab takes a fac- 
torizable form. 

Let us first assume that 



\z\) = aiki&i) 



\Z2) 



a 2 \a 2 b 2 



(27) 



where the local basis |aj) and \bi) are normalized. Ac- 
cording to relations {z\\zi) = {z 3 \z 3 ) = (zi\z 3 ) = 0, one 
can check that an arbitrary superposition of the two fac- 
torizable states \zi) and \z 3 ) shall still be a factorizable 
state of qubits A and B. So it can be concluded that 
the state \z 3 ) must take a form of \z 3 ) = a 3 |ai6 3 ) or 
\z 3 ) = a 3 |a 3 6i). Note there are the same relations for 
the states {\z 2 ), \z 3 )}. Combining these two restrictions 
gives that the state \z 3 ) must take the form 



W) = a 3 \aib 2 ) 



or 



*3 



a 3 |a 2 6i 



(28) 



(29) 



Correspondingly, the similar analysis is applicable to the 
states {\zt), \z 2 ), \z&)}. Knowing that the four states \zi) 
(i = 1, • • • , 4) are linearly independent, we can obtain the 
form of the state I24): 



\z 4 ) = a 4 |a2^i) 



or 



\Z4 



ot4\a 4 b 2 



(30) 



(31) 



corresponding to the two cases of (|28|) and (|29|), re- 
spectively. Finally, according to the relation (z\\z 2 ) = 



3 



(Z3I-S4) = 1, the complex coefficients of the states must 
satisfy laia^l = la^a^l for both of these two cases. Now 
it can be shown directly that the operator Fab takes a 
factorizable form of local filters of qubits A and B. 

Up to now, we have presented a distinct expression for 
the associated operator Fab an d shown that the local fil- 
ters /a and /b can be obtained through factorizing Fab 
straightforwardly. To complete the argument, we need 
to calculate the probability Pt that the LOCC action 
succeeds. Directly, 



Pf=tr(pfA®fl) = detfAdetf B J2( x i\ F AB\xi)- (32) 



It can be found, from Eqs. that the coefficient 

det /a det fs is equal to the minimum eigenvalue A^ in of 
the associated operator Fab ■ Thus we have 



P f = X^JrR(p). 



(33) 



Illustration. — We now apply the results derived above 
to treat the following four-parameters family of mixed 
states, 

p = p 1 \ip 1 )(tp 1 \ + P 2^2)(H +j>a|00)(00|+p4|ll)<H|. 

(34) 

Here 

= a|01) - 0\1O), = /8|01) + «|10) (35) 

and a, f3 are two positive real numbers satisfying a 2 + 
f3 2 = 1. The four states |^i), \tp 2 ), |00> and |11) com- 
pose the orthogonal decomposition of p and the corre- 
sponding probabilities are assumed to satisfy pi > P2 
and P3 > pi. Exploration of the states ([54]) should be of 
universal value: it comprises lower rank states as some of 
the coefficients pi reduce to zero, so the results acquired 
should predict the various consequences mentioned above 
for these reduced cases. We demonstrate these expecta- 
tions in detail below. 

A direct calculation gives the Wootters representation 
of the state p with 



1*2) 
1*3) 

\Xi) 



icos9y/pi~\ipi} 
sin0 1 /pi"|V'i) - 



- ism9^/p~2~\ip 2 ) 

■ CO8 y/p^\lp 2 ), 



V2 



(VPilOO) 



Pl|ll) 
Will), 



(36) 



and their tilde inner products are 



1/2 



Ai = [a 2 /3 2 {pi -p 2 ) 2 +P1P2 
\2 = [a 2 l3 2 { Pl -p 2 ) 2 +Pi P 2] 1/2 

A3 = A 4 = 

The parameter 9 in (|36|) is given by 



a(3{pi 
a/3(pi 



P2), 
P2), 



(37) 



= arctan 



VPlP2"(" 2 - P 2 ) 



a 2 (3 2 {pi -P2) 2 +PiP2} 1/2 + a(3[p\ +p 2 )' 

(38) 

Provided the presumption that the considered state p has 
nonzero entanglement of formation (thus nonzero concur- 
rence), the parameters a, /3, and piS should satisfy 



C(p) = 2[af3( Pl -P2) - VP^\ > 0. 



(39) 



Now according to the results derived in the preceding 
section, the local operation of the optimal transformation 
for the rank-4 state p can be uniquely determined. The 
factorized form of the associated operator Fab can be 
obtained by a step-by-step calculation from the proposed 
approach. The final form of it can be shown as 



1 



F AB = fc 2 (^J^|0)(0| + |l)(l|)®(|0)(0| 



1 



where 



(y/pi~a cos 9 — y/p2~f3 sin ( 



|1)<1|). 
(40) 



(41) 



The local filters j 'a and fs can be easily determined from 
E qs. (^) -(j^). For example, as the parameters satisfying 
yPifpi < yVz/pi < k 2 , they can be shown as 



/A = i(^) 1/4 |0)(0| + |l)(l|,/ s = 
k P3 



k Pi 

(42) 



and for the case k 2 < \Jpilp3 < \J P3 /p± , they are 

f A = |0)<0| + fc(^) 1 /4|i)(i|, /s = fc(^)l/4|0)( | + |1)(1|. 

Pi P3 

(43) 

The probability of the transformation is given by 

P f = ll { I" 2 / 3 ' tPi - V2? + PiP2] 1/2 + VPipl} (44) 



or 



P f = 2k 2 \[a 2 f3 2 {p 1 -p2) 2 +piP2] 1/2 + ^mpl) (45) 

corresponding to the above two cases, respectively. The 
transformed Bell diagonal state can be worked out from 
Eq. (|l|). Asp 2 = 0, it can be shown simply 



p' = -[2pia/3|¥_)<¥_| + ^1(1$+ )<<!>+ 1 + |$_)($. 



(46) 



where the normalized factor N = 2p\(x{5 + 2^fpip~i and 
I*-) and |<I>±) are Bell states given by 



4 



|*-> = ^(|01>-|10», 
|$±) = ^(|00)±|11)). 



(47) 



The state p' has the maximal extractable concurrence 



C(p') = 



Pia(3 - ^JpzPi 



C{P) 



trR(p) 



(48) 



P2=0 



which verifies the formula (12) presented in Ref. |12|. 

Now let us give a discussion for the degenerate cases. 

1. First consider the situation as the parameter p^ 
approaches to zero asymptotically and the limit of it cor- 
responds to a rank 3 density matrix with A3 = (or a 
rank 2 state with A 2 = as P2 = 0). In such a case, 
the relation \Jp±/ps < k 2 < \/pzjp± exists and the local 
filters Ja and Jb can be obtained from Eq. (^) as 



!a = -^) 1/4 |0)<0| + |1)<1|,/b = fc(^) 1/4 |0)( 
k p 3 p 3 



|1><1| 
(49) 



One can see in the limit of P4 — > 0, Ja and fs reduce 
to projection operators of pure states of the two qubits, 
respectively, and the probability of the manipulation 



Pi 




{[a 2 f3 2 ( Pl 



Pi) 2 +PiP 2 } 1/2 



(50) 



becomes zero. This confirms the assertion that any two- 
qubit density matrix with A„ = cannot be transformed 
into Bell diagonal state efficiently by invertible local ac- 
tions. 

2. Another notable situation is the rank 2 case with 
P3 = Pa = 0. In this case the expression given by ([4(i| ) 
appears to be not valid any more. The analysis proposed 
previously predicts that there exist various solutions for 
the local operator, and they can be obtained by estab- 
lishing the bases {\x 3 ), |x 4 ) } as 

\x 3 ) = *(ri|00) + r 2 |ll», la*) = n|00) - r 2 |ll), (51) 

where t\ and r 2 are complex parameters and the multi- 
plication of them satisfies tit 2 = |tit 2 | > 0. Similarly 
the factorizable form of the associated operator Fab can 
be worked out 



Fab = k 2 ( 



2, 1 N 



k 2 



|1)(1|)®(|0)(0| + ^M|1)<1|). 

(52) 



In comparison with Eq. (|40|), the difference is that the 
ratio |t 2 |/|ti| in ( [52] ) can be an arbitrary positive number, 
which implies the multiplicity of the local action. As the 
parameters are chosen that |ti| = |t 2 |, the local action 
has the optimal probability 



P f = 2[a 2 p 2 ( Pl -p 2 ) 2 +PiP2] 1/2 x min{fc 2 , i}. (Wi 



In conclusion, we have explored the optimal entangle- 
ment manipulation for a single copy of two-qubit states. 
By virtue of a particular representation of the two-qubit 
state, an explicit form of the LOCC action has been de- 
rived to transform the state into Bell diagonal state and 
the probability of the transformation was also given. The 
results display that the local operation of such a process 
is unique for the density matrices with the rank n — 3 
and 4, whereas it can be constructed variously for that 
with n = 2. 
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